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DIST M>u-riONS or rvwcrtow* o, «»h»om v**iabi.*s 
Thus independence is not required for the first result. Now 

ear - siri^*' " "tf) 



where 



(6-18) 



«.» SYNTHETIC SAMPLING (MONT6 CARLO TECHNIQUE) 

The methods discussed so far to find the distribution of 

Y = g(.Xi,Xi *".) 

have been approximation, or too involved to be pacticU for targe P^™' 
In this section a wry simple and intuitively satisfying nunhod is F^-^ 
Its only drawback* are that it require* a digital computer and fMU»*lF«»- 
metr* r«Sts are not obtained. thus limiting the W^ty te >V"*™*- 

It is assumed thai Y *-> is known and that the J™» 

YW * known- Now if a sample value of each random vanabte were 

Suld be computed [say * = |d., • - v^JJJ? * x 1 

sampleyalues were chosen for the random vanables (say *i - *n *- 

then y, - g<«„, *i. **0 «*uld be computed. 

If one had the time one could compute many such sample 
The computer actually supplies the speed that makes jnany s ^«^ 0 ™ 
JSS There is just one problem. How does the computer select the 

different values of X u Xt, ••• i X*1 j--,ii mt i ni , hh«A 0 

If each of the random variables had a uniform distribution between w 
and" " umbers for each random variable could be chosen ftom a ^ble of 
^dom numbers. Actually, computer routines generate pseudorandom 
cumbers which may be used. . x 

Consider the following case. Let the random var tables *i. ... , 
be independent and each uniformly distributed between zero and one, Jind 
kt Sr^JfA. ^beaknownfunction. Then thecomputer program 
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*.» S f t"*TM <T< C SAMPLING (MCNTA CARLO TECHNIQUE) 

to compute an approximation to the distribution of Y consists of the follow. 



g basic steps. 
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If a plot is desired it may be convenient to plot a usual bar chart or 
histogram. This method simply calls for breaking the range of Y into say 
30 mutually exclusive cells of the same size and plotting vertically the number 
of samples that fell into that cell. (See Chapter VIII for more details of 
histograms.) 

The case of uniformly distributed variables was considered. Now let X 4 
have a distribution function To obtain a random sample of X i the follow- 
ing procedure may be used. Select a random sample of V which is uniformly 
distributed between 0 and 1. Call this random sample i^. Then F^~ l ( u D' 1 ^ 
the random sample of X (see Example 6-7). 

For example* suppose that X\% uniformly distributed between 10 and 20. 
Then 



-1. 



x < 10, 

10 < x < 20, 



This is shown in Figure 6-29* 

Notice Fx^Hu) = lQu + 10, Thus if the value .250 were the random 
sample of U % then the corresponding random sample of A* would be 12.5. 
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As another example suppose that X is normally distributed with mean 0 
and variance 1. Then a random value of .250 for u J»uW ■ e «^* * 
_ 67 for a value of A'. This result follows from a table of the normal distriou- 
tioo function- Practically, most computers automatically generate atandard 
oorma ] random variables. If this is the case then random sample of ^ normal 
xandom variable Y with mean j/auda variance a- may be generated from a 
standard normal random variable X by recalling that 



X « 



r-i» 



or 



(6-19) 



Equation 6-19 can be easily checked by finding the characteristic fnnction 
of Y 

The only difference when the random var Ubles are dependent is that the 
dependence 'must be taken into account when the random sample* are ge ner- 
^ We assume that the dependence is expressed in terms of conditional 
distributions. If this is not the case the joint distributions can always be 
reduced to the required conditional distributions- 

For illustration consider three dependent random variables: *i» »™ 
X* We first generate a random sample of * by the same methods disco^ed 
above. Call this sample * u . We next generate a random ^P lc / f ^f^f 
^ IX _ t by the same method used before. Call this sample * 12 . Then we 
u£ x ^ to generate Thus nothing changes except that the 

. condiUon^tribulioa functions are used in generating the randomsamples. 

With a fast digital computer thousands of simukUons can be run xn 
reasonable limes. Monte Carlo solutions often involve 10,000 or more 
simulations. An example is given in the next chapter. 
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6-9 SUMMARY 

The purpose of this chapter was to consider the important engineering 
problem of finding ihc distribution of Y g(Xi* - - . , X n ) where the 
distribution of the JK/e is known. 

First the problem of Y *» g{X) was considered and this problem was 
solved by (6-2). Examples were given to illustrate its application. 

Next Y ™ -*« was considered, and it was shown that in ihc case 
of independent random variables the solution involved convolution of the 
density functions or multiplication of the characteristic functions. The case 
where Yit a linear combination of the X t f s was shown to be only a slight 
generalization of this problem. 

The general problem was then considered and although a general method 
of solution was outlined, the difficulty of solution was illustrated by an 
example* and approximations were suggested. 

Two approximations for the general solution were described. First a 
Taylor series approximation, moments, and the central limit theorem were 
used. Then a Monte Carlo method was suggested. 

References Bl, D3, D4> and P2 provide additional reading. 



6-10 PROBLEMS 

1. X has a normal density function with mean 1 and variance 2. 

I. 

Find the density of Y. 
2* Trie power F dissipated in a rooster is P = jfijL Assume A - 2 and / is a 
random variable with a normal density. 

1 -i* 

Ffcd the density function of/*. 
3* The output of a full wave rectifier is Y - \X\. Find the density function of Y 

when XXm a uniform density from -1 to -H- 
4. The output of a square law detector is 

y = or\ o > o. 

Find the density function of Y in terms of the density function fx ofX. 
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7-S INTRODUCTION TO TO*.**AWCB STUDIES 

is usually a good assumption* When the standard deviation of the pans axe 
known, then the mean and variance of the output can be computed using the 
approximation developed in the last chapter. 

Again something must be assumed to describe the output distribution and 
the probability of being out of tolerance. It is suggested that if enough 
variables are involved and the function is approximately a linear combination, 
then a normal density may he assumed. Then the probability of being outside 
tolerance limits can be computed* 

We now show an example of a tolerance problem whieh is solved using a 
Monte Carlo approach. 

Example 7-2 

This example, a simplification of a problem that actually occurred in 
practice, was worked by the method shown, and the results actually obtained 
in manufacture corresponded with the theoretical results. 

The simplified version is shown in Figure 7-4, The bar will fit within the 
bracket if Y 4 < Y> and there will be interference (or no fiO if Kb < >V 1« 
the actual case Y a and Y» involved 41 dimensions and the configuration was 
more complicated than simply a sum of lengths. Actually the equations for 
Y a and Y k involved arcs and angles, thus various trigonometric functions 
were involved. 

The problem was solved by finding the probabihry distribution 01 Z = 
n - Y A by Monte Carlo sampling. Note that if Z > 0 there is no problem, 
while if Z < 0 there will be no fit and the parts cannot be assembled. 

The problem arose because just as production and assembly were about 
to start it was discovered that interference was possible. (As customary, the 
intent was to tighten the tolerances on each part until no interference is 
possible at the worst case, but the designer made a mistake in his worst case 
calculations.) Then the question was, what is the probability of inter- 
ference? If it is low enough then it would be better to have a few that would 
not fit, rather than wail and spend the extra money to redesign some of the 
parts* 

To find the probability density of Z via a Monte Carlo technique* one 
must have Z =g&x *n) and must know the joint distribution of 
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^ The function * was found from the drawing* using trigonom- 

Jf 1> - - . ♦ The runcuon ^ w» i simply to illustrate the 

etry* A pact of the equation is shown m Figure 7-5 simpy 
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Figure 7*5 

forn, The various dimension, were turned to be ^^fJ^Stoa 
Ufccty between their wpper and lower tolerance WIS. The result or 

71 time, in 80<H> simulation to ? t ^JS^ ^^cTin^fcrencc 

well with the prediction that 71/S00O would not fit. 
Summary 

Tolerances of paramewrs eon be combined by ™»*^$£g£Z 
Taylor ser* -P^^^^^ 

tolerances w assign to what parts. 
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